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Abstract 

We show that a central l/r n singular potential (with n > 2) is renormalized 
by a one-parameter square-well counterterm; low-energy observables are made in- 
dependent of the square-well width by adjusting the square-well strength. We find 
a closed form expression for the renormalization-group evolution of the square-well 
counterterm. 
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1 Introduction 



The study of singular potentials in quantum mechanics is almost as old as quantum me- 
chanics itself jl). Physically, singular potentials pose problems because the force between 
two particles, represented by the potential, does not uniquely determine the scattering 
problem @. Here we will focus on singular potentials of l/r n type, with n > 2. Classi- 
cally, particles subject to such a force fall to the origin with an infinite velocity. In the 
quantum theory, the wavefunction oscillates indefinitely on the way to the origin, allow- 
ing no way of specifying a linear combination of solutions [[|. Of course, in any physical 
situation described by a singular potential, the potential is intended as a description of 
long-range behavior, so there is a sense in which the pathologies which occur near the 
origin are irrelevant to the physical problem. This should remind the reader of the infini- 
ties encountered in quantum field theory which are cured through renormalization. This 
analogy with field theory has provided an important motivation for the study of singular 
potentials ||. 

If a singular potential itself is not sufficient to determine the scattering problem, one 
might be tempted to classify the singular potentials as nonrenormalizable and abandon all 
hope. This point of view is now outdated. In the modern version of the renormalization 
paradigm a low-energy system with a clear-cut separation of scales can be described 
by an effective field theory (EFT) involving explicitly only the long-wavelength degrees 
of freedom, and organized as an expansion in powers of momenta [|[]. The short-range 
dynamics can always be treated as a set of local operators. In the present context, the l/r n 
potential represents the long-distance part of the potential. Local operators in momentum 
space correspond to delta-function interactions in coordinate space. The essential point of 
EFT is that the details of the short-distance physics are not of importance to low-energy 
scattering. Hence one can simulate the delta function in an infinite number of ways. The 
simplest choice of a "smeared out" delta function is a simple square well. With a singular 
potential representing a given long-distance force, and a square well representing unknown 
short-distance physics, an interesting question is whether one can obtain an EFT with 
well defined low-energy scattering observables, which are to a specified degree of accuracy 
insensitive to the short- distance physics encoded by the square well. It is the purpose 
of this paper to explore this issue. Note that we do not attempt to renormalize the 
coupling strength of the singular potential itself ||. In the physical problems of interest, 
the coupling strength is completely determined by the long-distance physics so there is 
no freedom to renormalize this parameter. 

By way of physical motivation we note that the singular potentials of l/r n type are 
of great current physical interest. The special case n = 2 is relevant to the three-body 
problem in nuclear physics @ || . This case is also relevant to point-dipole interactions 
in molecular physics M . The case n = 3 corresponds to the tensor force between nucleons 
and is at the heart of nuclear physics. The issue of the proper renormalization of this 
potential is an essential ingredient of the intense ongoing effort to develop a perturbative 
theory of nuclear interactions fll0|. The interaction between a charge and an induced 



dipole is of type n = 4 [11 1. The case n = 5 is a perturbative correction to the tensor 
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force in the nuclear potential QTOJ . Both n = 6 and n = 7 correspond to van der Waals 
forces, of London |T2| and Casimir-Polder [IBJ type, respectively. 

This paper is organized as follows. In section 2 we set up the quantum mechanical 
scattering problem of two particles subject to a l/r n potential with a square well. In 
section 3 we consider the marginal n = 2 case in some detail. The pure singular potentials 
n > 3 are considered at zero energy in section 4. In section 5 we make use of the 
WKB approximation to generalize our results to non-zero energy and to estimate the 
errors associated with the renormalization procedure. We discuss the applicability of a 
perturbative expansion for singular potentials in section 6. Our numerical analysis, for 
the case n = 4, is discussed in section 7. We discuss and conclude in section 8. 



2 The l/r n potential with a Square Well 

We consider two particles of reduced mass M interacting in the S wave with a singular 
potential that goes as l/r n , n > 2. This potential has a scale that sets its curvature, 
ro] this is the characteristic scale of the long-distance physics. The strength of the long- 
distance potential is governed by a parameter \^j2Mr\. To obtain well-defined solutions, 
we need to regulate the potential by introducing a cutoff procedure. Since we are posing 
the problem in coordinate space, we do this through a cutoff radius R = 7Zr , 1Z<1. 
We expect that the solutions will depend sensitively on 71, that is, on the short-range 
physics. We simulate a short-range delta-function interaction by a square well of this 
radius and with depth Xs/2Mtq. The problem will be correctly renormalized once we 
are able to vary 1Z (inasmuch as 1Z<1) and simultaneously A$ in such a way as to 
keep observables (say phase shifts) invariant. The corresponding constraint \$ = ^s(7Z) 
represents the renormalization-group flow of the contact interaction []. We will see that 
this short- distance physics is represented by the running coupling 

H n (K) = ^\ S (K)K. (1) 

We thus take as our potential 

where f(x) is a regular function of x near the origin with /(0) = 1, /(l) = 0(1). Notice 
that Ag, \l > correspond to purely attractive potentials. In terms of x = r/r , the 
Schrodinger equation for the wavefunction u(r)/r at an energy E = k 2 /2M = rj 2 j2Mr\ 
is 

u"{x) + {rj 2 + X s )u(x) = x < K 

u"(x) + (ri 2 + \ L ^)u(x) = x>1Z. { ) 

7 Here we mean that there is a "group" of transformations on the cutoff 1Z which leaves observables 
invariant. 
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We will consider the simplified case f(x) = 1 until section 5. 

There is a very simple argument for classifying singular potentials which we will repeat 



here [fL^] . In the vicinity of the origin the uncertainty principle dictates that the kinetic 
energy scales like x~ 2 . Therefore in a system described by an attractive singular potential 
alone, the Hamiltonian of the system is given by the sum of the kinetic energy and the 
potential energy, — \ix~ n . Note that for a Coulomb potential, n = 1, and for a sufficiently 
weak n = 2 potential, the Hamiltonian is bounded and therefore the Schrodinger equation 
has a unique regular solution. Clearly for n = 2 and Xl sufficiently strong, the Hamiltonian 
is unbounded from below. Furthermore, when n > 3 the Hamiltonian is always unbounded 
from below. Hence, an attractive singular potential alone is meaningless in the vicinity of 
the origin; the unboundedness of the Hamiltonian represents the onset of short- distance 
physics whose effect must be included in the potential. 



3 The Marginal Case: n = 2 
3.1 The k = Solution 

Consider first the zero-energy solution u(x; 0). The Schrodinger equation for x > 1Z is 

u"(x;0) + ^u(x;0) = (4) 
x l 

and the general solution is 

u(x; 0) = Ax^ +1 + Bx^~~< (5) 

where 7 = \Jl/A — Xl- For A^ < 1/4 the solution is well known |12| and will not be 
further considered in this paper. On the other hand, for A^ > 1/4 we can define 7 = iu, 
and the general solution is 

u(x; 0) = v^ccos (z^log (x) + fa) (6) 

where v = \JXl — 1/4, fa = (\ogA/B)/2i and we ignore the overall normalization. Both 
of the linearly independent solutions of Eq. ([]) vanish as x — > 0, and oscillate indefinitely 
on the way there. There is no obvious way to determine a unique linear combination 
of solutions; i.e. fix fa. This is the fundamental problem with singular potentials in 
quantum mechanics. Renormalization theory tells us that this sickness is to be expected 
and arises from probing arbitrarily short-distance scales, where the true potential no 
longer has the form 1/x 2 . The cure is to cut off the long-distance potential at a radius 
1Z and introduce a simple parametrization of the unknown short-distance physics, since 
low-energy observables cannot distinguish between the schematic, parametrized potential 
and the true potential at short distances. We choose a square well for simplicity, but we 
emphasize that any choice of function is equally valid. 
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3.2 Matching to the Square Well 

The solution in the interior region, x < TZ, is straightforward. It is sufficient to consider 
an attractive square-well potential, Xs > 0. Matching logarithmic derivatives at the 
boundary x — TZ gives 

\ s coty/%H = - j--i/tan(z/log(ft) + 02)}. (7) 

If we vary TZ and Xs(TZ) as given here, the zero-energy phase fa will not be affected. Eq. (0) 
is transcendental and therefore rather cumbersome. However there are two regimes where 
an analytical expression can be found. The first is when cot \f\~sTZ is large. This is the 
generic situation as TZ — > since the the right-hand side of this equation blows up, except 
where z/tan(^log (TZ) + fa) = 0. We can then find an approximate solution by writing 
y/X~sTZ = ran + e where m is an integer and e is a small number. Inserting this into Eq. (|7]) 
and keeping the leading order in e we find 



H 2 (TZ) = mir 



sin (u log (TZ) + fa — arctan j- 



sin {v log (TZ) + fa + arctan 



Close to the zeroes of the right-hand side of this equation we can write instead \f\~sTZ = 
(m + 1/2)tt + e and following the same procedure we find, to leading order in e, 

H 2 (TZ) = (m + 1/2)tt - + * [\ - ^tan (ulog (TZ) + 2 )} . (9) 

A numerical solution of Eq. (0) shows that the approximate solutions Eqs. are very 
good within their ranges of validity. The expression in Eq. @ interpolates between two 
successive branches of Eq. (g) (see Fig. [I]). 

The three-body problem with short-range interactions is known to be equivalent, in 
the ultraviolet regime, to the 1/r 2 potential 0. The role of the interparticle distance 
r is played in the three-body problem by a collective coordinate that vanishes when 
the three particles occupy the same point in space, and the analogue of the three-body 
force is the short- distance potential. Since renormalization depends only on the short- 
distance behavior of the theory, it is not surprising that the renormalization of the three- 
body problem requires the presence of a three-body force (J7|. Using momentum-cutoff 
regularization the running of the three-body force was found to follow Eq. (H) very closely, 
even where this formula predicts a pole. Evidently, for values of TZ where Eq. ([]) takes 
over in the 1/r 2 problem, the three-body force continues to follow Eq. (§) and seems to 
reach arbitrarily high values. Also, no evidence of multiple branches was found in the 
three-body problem. These discrepancies between the three-body and the 1/r 2 problem 
may be due to the fact that not all aspects of the renomalization-group flow are universal, 
in the sense of being independent of the particular regulator used. 
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Figure 1: The running coupling for the n = 2 singular potential. The solid lines are given 
by Eq. (Qj and the dashed lines are given by Eq. ^). The bold lines are a numerical 
solution of Eq. fffy. 



3.3 The Full Solution 

In the case n = 2, the Schrodinger equation can be solved exactly for all energies. The 
solution is 

u(x; rj) = \fx [exp (ia) J iu (rjx) + exp (— ia) J_ iu (r}x)] (10) 

where the J±i U are Bessel functions, and a is to be fixed by a boundary condition. For 
small x we find 

u(x; rj) = \[x cos (flog (xrj/2) + a — Imlogr(l + iu)). (11) 
Matching to Eq. (^) gives 

2 = a + z/log77/2-Imlogr(l + iz/). (12) 

Since 2 is, by construction, energy independent, a is energy dependent. 

We can now look for solutions with rj = in which fall off exponentially at large x. It 
follows from Eq. ( |l0l) that 

1 7T VK 

u(x; k) — > - exp (i— ) cos(a + *~^~) ex P { KX ) + C exp (— kx) (13) 

where C is an energy- dependent coefficient. The bound-state solutions then correspond 
to a(r]) = (m + l/2)7r — iz/7r/2, with m an integer. Comparing with Eq. (|12"D gives the 
bound-state spectrum 
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E m =-^ e X p 2 fe + tologr(l + H-(m + l/2M (M) 

Once 02 is fixed by a single bound-state energy, all other energies are predicted 0. 
Adjacent bound-state energies are related by 

/ 71 \ /ir\ 

= exp( ). (15) 

Hence we see that the periodicity in the running coupling H2(TZ) is associated with the 
accumulation or dissipation of bound states near the origin. 

One can also fix <f) 2 to a scattering observable, like the scattering length or the phase 
shift at a given energy. Unfortunately, as for the Coulomb potential, the n = 2, 3 singular 
potentials suffer infrared problems at low energies, and therefore scattering lengths can 
be defined only if an infrared cutoff is imposed ffl. 



4 Pure Singular Potentials: n > 3 
4.1 The k = Solution 

The exact zero-energy solution for n > 3 is well known ||. Defining z = \fXlx l ~ n l 2 /|1 — 
n/2\ and (f)(z) = u(x; 0)/y/x, for x > 1Z Eq. (^|) becomes an ordinary Bessel equation: 

f(z) + U'{z) {n _\ )2z2 ) <!>{*) = 0. (16) 

The solution is 



u(x; 0) = \fx 



A n Jl/(n-2) I — 



^—x l ' n / 2 )+B n J. 



n/2 



n J -l/(n-2) 



1 — n/2 



—x x ~ nl2 



which is a linear combination of Bessel functions. For small x we can write 

^7 



u(x; 0) = x n/4 cos 



n/2 - 1 

where we have set the constant pref actor to unity and 
, w 1 ( B n 

^ = -T(^2) +dog ( 1 + a: 6XP 



1 + C(x 



n/2-l> 



Z7T 



(n-2). 



This solution exhibits precisely the same pathologies as Eq. 



(17) 



(19) 



3 In the case n = 4, the Bessel functions are of half-integral order and Eq. (|18| ) is exact for all x. 
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4.2 Matching to the Square Well 



We proceed as in the case n = 2. Again we have a square well in the interior region. 
Matching logarithmic derivatives at the boundary x — 1Z gives 



A s cot J\~ s n =^--(hA 112 tan ( -^^n 1 ^' 2 + n ) (20) 
v ATI \n n j \n/2-\ Y J v ' 

where we have neglected 0{lZ n ^ 2 ^ 1 ) corrections to the wavefunction at x > 1Z. The phase 
4> n is physical and can be traded for the scattering length (for n > 3), as will be seen 
below. If we vary 1Z and Xs(TZ) as given here, the phase n will not be affected. We 
proceed as we did before in the n = 2 case and find, in the regions where the right-hand 
side of Eq. (pD|) is large 

HJK) = run J 1 j—= r \ . (21) 

\ 1 - n/4 + ^T L n l ' n / 2 tan (^ll 1 ^/ 2 + ' ' ' 

In the other regime, where the right-hand side is close to a zero, we have 

H n {U) = (m + J)tt - . I r- fn/4 - ^n 1 ^ 2 tan f^W-™/ 2 + 0^ V (22) 
2 (m + i)7r \ v \n - 2 // 



A numerical solution of Eq. (p0|) in the case n = 4 shows that the approximate solutions 
Eqs. (|2l|j22] ) are very good within their ranges of validity. The expression in Eq. (|2"2"[) 
interpolates between two successive branches of Eq. ( pip (see Fig. |^) . 



The scattering length can be found from the zero-energy wavefunction for n > 4 pA 
For instance, we find the n = 4 scattering length 



04 = r y Aitan04. (23) 
It is evident that 04 determines the phase 04. 



5 The WKB Approximation 

There is an important shortcoming in what we have done so far. Defining zero-energy 
scattering is not sufficient to guarantee correct renormalization. We want physics at all 
energies 77 < 1 to be cutoff independent. It is clear that the above procedure could in 
principle be repeated at each energy by allowing energy dependence in A5. Fixing <ft n at 
one energy and then predicting other energies will only result if the scale of this energy 
dependence is much slower than l/2Mr 2 , so that, to some accuracy, As can be taken 
to be energy independent. Otherwise, an infinite number of parameters (the strength of 
arbitrarily-many-derivative contact interactions) would have to be known in order to have 
predictive power. This shortcoming can be removed using the WKB approximation. We 
can also consider the more general case, f{x) ^ 1. 
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Figure 2: The running coupling for the n = 4 singular potential. The solid line is given 
by Eq. and the dashed line is given by Eq. (WW). The triangles and stars are extracted 
from numerically solving the Schrddinger equation and coincide exactly with the numerical 
solution of Eq. 



5.1 The WKB Criteria 

We now keep f{x) arbitrary and consider the region x > 1Z in the limit x — > 0. A 
particularly well-suited approximation in this limit is the WKB approximation, which is 
valid when the wavelength A is small compared to the characteristic distance over which 
the potential varies appreciably. That is 



^l£l = l£p"(S-V(r))]-V»|<l, 



(24) 



which translates into the constraint 



n 



-J\ L f(x) > x 



n/2-1 



(25) 



in the small-x region^. Clearly this condition is satisfied for all n > 2 as x — > 0. In 
the marginal case n = 2 this condition is satisfied only for a sufficiently strong potential. 
Therefore the WKB criterion parallels the general argument given above based on the 



boundedness of the Hamiltonian. The general WKB solution El is 



u(x; rj) 



> + WO*) 



-1/4 



cos I J dx I if + 



1/2N 



rr> 1 1> i If 1 ryl ■ 

Ju I \ J XQ \ Ju 

where Xq is a constant of integration. For l/>£ this reduces to 



(26) 



The WKB approximation is also valid at large x and finite rj provided that 77 nX^x n jl 
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u{x; 0) = x n/4 f~ 1/A (x) cos (\/^ jf dx'x'- n/2 f 1/2 (x')^ . (27) 

In the limit 1Z < x -C 1, we can set f(x) = 1 (keep the leading term in a power series in 
x). We then recover, for n > 2, 

0) = x™/ 4 cos (^^Y^ 1 "" 72 + 4>n\ (28) 

where n = — y/\LXo~ n ^ 2+1 / (I — n/2). The case n = 2 is also recovered if one takes 
— > \l — 1/4. Therefore, we expect our conclusions about the renormalization of the 
singular potentials to be valid for the more general case f(x) ^ 1. 

5.2 The Leading Energy Dependence 

We now show that the zero-energy solution is in fact sufficient to remove cutoff dependence 
at all other low energies. The crucial point is that, in the intermediate region TZ < x <C 1, 
for the energies of interest, the potential energy is much larger than the total energy, and 
we recover the zero-energy case. This can be made more precise using WKB again 0. 
We write the wavef unction for any x > 1Z as 

u(x; rj) = A(x; r))u(x; 0). (29) 

Then A(x; if) obeys 

d 2 A(x) din u(x: 0) dA(x) 9 . . . . , 

—\± + 2 ^ v ' ; } ! + V 2 A(x) = 30 

dx z dx dx 

which depends only on the zero energy wavefunction. Now, since for TZ < x <^ 1, 



dlnu(x; 0) 



dx 



» 1, (31) 



A(x) can be written 



A(x)=A m {x)+A {1) (x) + ..., (32) 



where 



dx ' dx 2m'(i;0)' (0) ' 



We then find the leading energy corrections 



^HAo,{l4/>'^ + ..] (34) 
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We see that, in the intermediate region, the energy dependence of the wavefunction (|29|) 
is determined by the zero-energy wavefunction u(x;0). If the phase of u(x; 0) has been 
fixed, the phase of u(x; rj) is fixed, and scattering observables can be predicted at low 
energies. 

5.3 Error Estimates 

The fact remains that our arguments are all at short distances where the WKB approxi- 
mation is valid. This is, of course, the opposite of the EFT limit which interests us. One 
may wonder whether cutoff effects can be amplified when propagating the wavefunction 
from short to long distances. We will see now that this cannot occur and in turn find 
an estimate of the cutoff error associated with the scattering phase shift. Usually, in 
perturbative EFT, the error is a power law in 1Z. Here we will find a more complicated 
functional dependence. 



By adjusting H n (7Z) as in Eqs. (£L],|3) we guarantee that two zero-energy solutions 
uk(x;0), utz'(x;0) corresponding to two different cutoffs 1Z, 1Z' <C 1 ~ 1 / 77 are identical. 
At finite values of 77, solutions obtained with different cutoffs will no longer be equal, 
but their difference can be easily estimated. Taking 1Z' < 1Z, the Schrodinger equations 
satisfied by u-ji(x; if) and u%i(x\ 77) are the same in the x > 1Z region so their Wronskian, 

W[u n ,u n ,}(x;ri) = u n (x; rj)u' n ,{x\ 77) - u' n {x\ rj)uw(x; rf) (35) 

is independent of x. At large distances (r ^> (A^/A; 2 )^/" - )), where the solutions are plane 
waves, W[u-ji,utii] is related to the phase shifts 8k,6-ji> obtained with the cutoffs 1Z and 
ft' by 

W[un,UK>}(r » (X L /k 2 f /n) ;v) = A n A n , V sm(5 n - 5 W ), (36) 

where An, An> are the amplitudes at large distances. These prefactors are easily estimated 
from the general WKB solution, Eq. (|26|) , in the region r ^> (A^/A; 2 )^ 1//n - ) , where the 
WKB solution maps to the asymptotic plane- wave solution (see Footnote 9). We find 
A n ,A n , -rr 1 / 2 . 

On the other hand, at the cutoff distance x — 1Z, W[uk,uii>] is estimated using our 
WKB formula, Eq. fl3|). We find 



where 



2 

W[u n , u n/ ] {TZ; 77) = W[u n , u w ] {1Z; 0) - ^-£(1Z; 0) + . . . (37) 



W[uTz,un')(TZ;0)=UTz(TZ; 0)^(ft; 0) 



u 



K 



,(K; 0) 



ti- 



ll 



(71; 0) 



o iz n/2 - 



-1 U K 



(3f 



and 
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£(1Z;0) = utz(K;0)utz>(K;0) 



uk>(K;0) _ u n (1Z;0) L n/2 -! %g) \ 
u' n ,(1Z;0) u' n (K;0) [ u' n (K)) 



+ W[ua,UR>](K;0) dx' 
Jo 



n 



u^(x';0) u n (x';0) 
u' n ,(x>;0) u' n (x';0) 



O 



x 




(39) 



We have included the error due to keeping only the leading zero-energy wavefunction 
in Eq. (|i~8|). Recall that we choose our fitting procedure to be energy independent, for 
example, by comparing the zero-energy wavefunction to the scattering length. It then 
follows that W[un,Ufi'](7l;0) = 0, by construction, for the full wavefunction, and from 
Eq. (KS) we have 



u> n ,(1l;0) u' n (1Z;0) 



O K n/2 - 



u n >(K;Q) u n (1Z;0) 
Using these constraints it is straightforward to find 



£(1Z; 0) = O [IT 12 



(40) 



(41) 



u' n {K) )■ 

If we assume that all oscillating functions of 1Z are of order unity for values of 1Z at which 
we fit observables, then £(Jl;Q) = OiTZ 311 ^ 2 ^ 1 ), which is small for all n > 2. Matching 
the Wronskians at large (r 3> (Ai/Zc 2 )^ 1 /" - ^) and short (a; = 1Z) distances then yields an 
estimate for the error in the phase shift: 



(42) 



where £(11; 0) is a function of 1Z whose complicated parametric cutoff dependence is given 
by Eq. (fil]). This shows that the renormalization procedure described here produces cutoff 
independent phase shifts, accurate up to order rf £(1Z; 0). 



6 The Weak Coupling Limit 

It might seem odd that the explicit dependence on the coupling constant is nonanalytic 
in the formula for the n = 4 scattering length, Eq. (|23|). Naively it would appear that 
nonperturbative effects are important at arbitrarily weak coupling. 

However, we know that this cannot be the case, since for weak coupling the scattering 
length should go smoothly to its square-well value. We would expect a perturbative 
description in the singular potential to be valid when the potential energy, — X L x~ n , is 
much smaller than the kinetic energy, x~ 2 . This leads to the condition 

r»r Af , or k^r^X^. (43) 
In effect, taking 4 from Eq. ( p0| ) we find, for Xl/1Z 2 1, 
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Figure 3: The exact zero-energy wavefunction (solid line) at small r compared with the 
wavefunction obtained in perturbation theory to leading order (small dashes), next-to- 
leading order (medium dashes), and next-to-next-to-leading order (large dashes). 



Leading order reproduces the square-well scattering length and the corrections are analytic 
in Xl- Hence, there is, in fact, no nonanalyticity near zero coupling in the presence of the 
square well. 

Of course, if the cutoff 1Z is taken at values where the oscillatory behavior of the 

i 

wavefunction has set in, 1Z < A£~ 2 , then there is no sense in which perturbation theory in 
Xl can capture the true behavior of the wavefunction. This is made clear in Fig. |3| where 
several orders in a perturbative expansion of the n = 4 singular-potential wavefunction 
are plotted against the exact singular-potential wavefunction in the short-distance region. 

7 Numerics 

In this section we analyze the 1/r 4 potential numerically. For simplicity, we take 2M = 
ro = 1; therefore, x = r and rj — k — y/2E. The long-distance potential is then completely 
determined by Xl which we take to be unity. We consider the "natural case", which 
is characterized by 0,4 ~ (Al) 1//2 , and the "unnatural case", which is characterized by 
a 4 » {Xl) 1 ' 2 . 

In Fig. [|we show phase shifts 8(k) in a natural case (5(0.1) = 0.1 and </>4 = —101.298) 
for various cutoffs. We see that, as anticipated, the low-energy phases are to a good 
approximation cutoff independent; cutoff dependence becomes more pronounced as the 
cutoff radius and the energy are increased. In the same figure we also plot the error 
analysis: the (log of the) errors |A5(fe)| = 5r — 5r> as a function of (the log of the) energy 
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Figure 4: Left: phase shifts 5 vs. energy in the case of a natural scattering length for n = 4. 
Various cutoffs are given with the square well tuned to give the same scattering length. 
The cutoffs are R = .01 (dots), R = .02 (triangles), R = .04 (squares), R = .08 (stars), 
R = .16 (diamonds). Right: logarithm of the errors \AS(k)\ = 5r — 8ri as a function of the 
logarithm of the energy for n = 4. The pairs of cutoffs (R, R') are: (.16, .08) (diamonds), 
(.08, .04,) (stars), (.04, .02, ) (squares), (.02, .01,) (triangles). 



for various pairs of cutoffs. We find that the errors scale as k 2 , as expected on the basis of 
Eq. ([42]) . In Fig. ^ we show the corresponding results in an "unnatural" case (5(0.1) = 7r/3 
and 04 = —98.954). Again we find that the errors scale as k 2 . 



8 Conclusion 

We have reconsidered singular potentials of the form l/r n with n > 2 from the viewpoint of 
modern renormalization theory. We have shown that the well-known pathologies near the 
origin are cured by a square- well counterterm which represents the effect of unknown short- 
distance physics. The renormalizat ion-group evolution of this counterterm has periodic 
behavior. The counterterm is not determined uniquely at any given cutoff due to the 
infinite number of branches of the renormalization-group flow, and one is allowed to 
freely jump from branch to branch at will without causing any change in the low-energy 
phase shifts (up to 0{k 2 ) corrections). 

The dependence of the number of bound states on the choice of branch is complex. 
Arguments similar to the one leading to Eq. ([14]) are valid for a generic singular potential 
as long as (i) the binding energy is much smaller that Xs/i^Mr^) (in order that the binding 
energy can be disregarded in the left-hand side of Eq. (p0|)) and (ii) the binding energy 
is much larger than Xl/ \2Mrl) (so the fact that f(r/r ) ^ 1 is unconsequential). The 



resulting spectrum shows a power-law distribution that is given by the WKB estimate fl"4" |. 
We see now that as 1Z —* and Xs(7V) — > oo the region of validity of the calculation 
sketched above grows and more and more bound states are created. If, in addition to 
keeping the phase shifts cutoff independent, one also demands that the number of bound 
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Figure 5: Same as Fig. [| for an unnatural scattering length. 



states also be fixed, the value of the counterterm \$ jumps down a branch at every cycle. 
We are then left with a periodic, limit-cycle behavior for the renormalization-group flow. 
This is a unique situation in field theory/critical phenomena, where the standard behavior 
corresponds to counterterms approaching either zero (asymptotic freedom) or infinity as 
the momentum cutoff goes to infinity |§. 

Naively it appears that the singular potentials have a nonanalytic dependence on the 
coupling parameter even at weak coupling. This would negate a perturbative description 
at weak coupling, a conclusion which must be incorrect. One might imagine some as 
yet experimentally invisible light particle which interacts at long distances via a singular 
potential (e.g. an axion). If the behavior of the wavefunction were such that there is a 
branch point at the origin of the coupling-constant plane, then nonperturbative effects 
would persist even for couplings of gravitational strength. We have seen that this nonan- 
alyticity is an artifact which is removed by short-distance physics encoded by the square 
well. 

Renormalization renders low-energy phase shifts cutoff independent up to 0(k 2 ) cor- 
rections. The cutoff dependence of these errors is not generally a power law as one expects 
in Wilsonian EFT. The renormalization-group flow introduces complicated oscillatory be- 
havior in the corrections, which nonetheless is small for judiciously chosen cutoffs. Our 
theoretical expectations of the error have been confirmed numerically. We expect that the 
methods developed in this paper will prove useful to those interested in the cornucopia of 
physical systems whose long-distance behavior is governed by singular potentials. 
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